COLLAPSING ESTIMATES AND THE RIGOROUS DERIVATION 
OF THE 2D CUBIC NONLINEAR SCHRODINGER EQUATION 
WITH ANISOTROPIC SWITCHABLE QUADRATIC TRAPS 



XUWEN CHEN 



Abstract. We consider the 2d and 3d many body Schrodinger equations in 
the presence of anisotropic switchable quadratic traps. We extend and im- 
prove the collapsing estimates in Klainerman-Machedon [25] and Kirkpatrick- 
Schlcin-Staffilani [23]. Together with an anisotropic version of the general- 
ized lens transform in Carles [3], we derive rigorously the cubic NLS with 
anisotropic switchable quadratic traps in 2d through a modified Elgart-Erdos- 
Schlcin-Yau procedure. For the 3d case, we establish the uniqueness of the 
corresponding Gross-Pitaevskii hierarchy without the assumption of factorizcd 
initial data. 



1. Introduction 

Bose-Einstein condensation (BEC) is the phenomenon that particles of integer 
spin ("Bosons") occupy a macroscopic quantum state. The first experimental ob- 
servation of BEC in an interacting atomic gas occurred in 1995 [1, 10]. Many similar 
experiments were performed later [9, 22, 29]. In these laboratory experiments, the 
particles are initially confined by traps, e.g., the magnetic fields in [1, 10], then the 
traps are switched in order to enable observation. To be more precise about the 
word "switch": in [1, 10] the trap is removed, in [29] the initial magnetic trap is 
switched to an optical trap, in [9] the trap is turned off in 2 spatial directions to 
generate a 2d Bose gas. The dynamic during the period when the trap is shifting is 
sophisticated. To model the evolution in this process, we use a quadratic potential 
multiplied by a switch function in each spatial direction for analysis in this paper. 
This simplified yet reasonably general model is expected to capture the salient fea- 
tures of the actual traps: on the one hand the quadratic potential varies slowly and 
tends to oo as \x\ — > oo; on the other hand, the switch functions describe the space- 
time anisotropic properties of the confining potential. In the physics literature, 
Lieb, Seiringer and Yngvason remarked in [26] that the confining potential is typ- 
ically ~ \x\ in the available experiments. Mathematically speaking, the strongest 
trap we can deal with in the usual regularity setting of NLS is the quadratic trap 
since the work [30] by Yajima and Zhang points out that the ordinary Strichartz 
estimates start to fail as the trap exceeds quadratic. 

Motivated by the above considerations, we aim to investigate the evolution of 
a many-body Boson system during the alteration of the trap. The N-body wave 
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function tp^(T 7 yN) solves the many body Schrodingcr equation with anisotropic 
switchable quadratic traps: 



id T ip N 

V'jvCo.yjv) 



\h^{t)iP n + 1 N n ?V(N? (y, - y,))V^ (1-1) 



2 

N 



n^o(yj). 

3=1 



3,1 



3,1 



(1.2) 



where rel, yjv = (yi,y2, — >yjv) <= R rlJV , V is the interaction between particles, 
and 

N N / „ / 2 \\ 

3 = 1 3=1 \l=l \ 

with the switch functions ?7/(t), Z = l,...,n. Throughout this paper, we only con- 
sider n = 2 or 3 and we assume the switch functions r/ l e C 1 (M ( | — > Mq ) satisfy the 
following conditions. 

Condition 1. 17; (0) = i.e. The trap is not at a switching stage initially. 

Condition 2. f) l is supported in [0, To] and Xo\/sup T |t?j(t)| < §. 

When the trap is fully on, Lieb, Seiringer, Solovej and Yngvason showed that 
the ground state of the Hamiltonian exhibits complete BEC in [27], provided that 
the trapping potential Vt rap (x) satisfies mi\ x \ >R Vtrap{x) —¥ 00 for R — > 00 and the 
interaction potential is spherically symmetric. To be more precise, let tp N0 be the 
ground state, then 



'JV.O 



\4>GP) 



IJ GP\ 



as N — > 00 



where 7^0 ^ s ^ ne corresponding one particle marginal density defined via formula 
1.3 and 4> GP minimizes the Gross-Pitaevskii energy functional 

( |V^| 2 + V trap (x) |^| 2 + 4™ |0| 4 )dx. 

Because we are now considering the evolution while the trap is changing, we start 
with a BEC state / factorized state in equation 1.1. 

However, tp N does not remain a product of one-particle states i.e. 

N 

^M^yw) ^ W^T.yj)^ > 

3=1 

for some one particle state <j>. Moreover it is unrealistic to solve the iV-body equation 
1.1 for large N. Thence, to observe BEC, we have to show mathematically that tp N 

N 

is very close to J^J 4>{t, yj), the mean field approximation, in an appropriate sense. 

3=1 

Notice that when <f> ^ </>' 

2 

N N 

n^ T 'yj)-n^( T '*>') - 2 " : " v :> x - 

3=1 3=1 
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i.e. our desired limit (the BEC state) is not stable against small perturbations. One 
way to circumvent this difficulty is to use the concept of the k-particlc marginal 
density 7^ associated with ip N defined as 



v P — c v 

l ( N(T,yi;y k ) = / ip N (r, y£, y N -l)^ N ( T ^ Yk, yN-l))dy N -l, yt,y k 



7>7lk 



(1.3) 

Another way is to add a second order correction to the mean field approximation. 
See [8, 20, 21]. 

In this paper, we take the marginal density approach and establish the following 
theorem. 

Theorem 1. Consider the 2d case when /3 G (0, |). Assume the interaction po- 
tential V is nonnegative and belongs to L 1 f] W 2,oc and the switch functions r/ l 
satisfy Conditions 1 and 2. Moreover, suppose the initial data has bounded energy 
per particle 



sup — {ip N ,H N (T)i> N ) 

N 



< CO. 



T=0 



where the Hamiltonian H n (t) is 



j=i \;=i \ y i-i J / i< 3 

If {7^} are the marginal densities associated with ip N , the solution of the N-body 
Schrddinger equation 1.1, and (j> solves the 2d Gross- Pitaevskii equation: 



id T 



0(0, y) 



0o (y). 



where H y {r) is the operator inside formula 1.2 and b n = J V(x)dx, then Vr G [0, T ] 
and k ^ 1, we have the convergence: 



7lv fc) (T,yt;yt)-n0(r, y ,)^(r, y ;.) 



as N -> 00. 



Example 1. We give a simple example to explain the switching process we are 
considering here: say 

r?i( T ) = C\ when t G ( — 00, ^], C2 when t G [1, 00), 

1 3 

77 2 (r) = C3 w/ien r G (—00, — ], C4 w/ien r G [—,00). 

Then our switching process contains the cases: turning off / on: C2 — / C\ = 
and tuning up / down: C\ < C2 / C2 ^ Ci. ,4s Zong as r/i(r) G C 1 ana" satisfies 
Condition 2, i] 1 can behave as one likes inside [i, 1]. Same comment applies to r\ 2 
too. Furthermore, Theorem 1 addresses the time intervals (— 00, 0] and [§,00) as 
weZZ, 5mce the equation is time translation invariant in these two intervals, we can 
use Theorem 1 separately in each sufficiently small time intervals. 
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Remark 1. Technically, one should interpret Conditions 1 and 2 in the following 
way. Due to Condition 1, we have a C 1 even extension ofrj l i.e. we define t]i(t) = 
??;(— t) for t < 0. The fast switching condition 2 in fact ensures that f3 l defined via 
equation 4-1 is non-zero in [0,T ] which is crucial in this paper. See Claim 1 for 
the proof. 

Remark 2. We assume [3 £ (0, §) to match Kirkpatrick-Schlein-Staffilani [23] in 
which the authors studied the n l = case. (3 = will yield a Hartree equation 
instead of the cubic NLS. 

(k) 

The approach with j y N ' has been proven to be successful in the rj l = and 
n = 3 case, which corresponds to the evolution after the removal of the traps, in 
the fundamental papers [11, 12, 13, 14, 15, 16, 17] by Elgart, Erdos, Schlcin, and 
Yau. Their program, outlined by Spohn [28], consists of two principal parts: on 

the one hand, they prove that an appropriate limit of the sequence yy N j as 
N — > 00 solves the Gross-Pitaevskii hierarchy 

(id t + \ - 7 « = b J2 B jlk+ i ( 7 (fc+1) ) ,k = l, n, ... (1.4) 

\ / j=i 

where Bjj~ + i are in formula 1.7; on the other hand, they show that hierarchy 1.4 
has a unique solution which is therefore a completely factored state. However, the 
uniqueness theory for hierarchy 1.4 is surprisingly delicate due to the fact that 
it is a system of infinitely many coupled equations over an unbounded number of 
variables. In [25], by assuming a space-time bound, Klainerman and Machedon gave 
another proof of the uniqueness in [14] through a collapsing estimate and a board 
game argument. We call the space-time estimates of the solution of Schrodinger 
equations restricted to a subspace of R" " collapsing estimates" . We can interpret 
them as local smoothing estimates for which integrating in time results in a gain 
of one hidden derivative in the sense of the trace theorem. To be specific, the 
collapsing estimate of [25] reads: Suppose solves 

« + i/w- 5^) 

there is C > 0, independent of j, k or u(' c + 1 )(0, xT/Tjf; x' fe+1 ) s.t. 



^ C 



fe+1 

n K v 



(1.5) 



\3 = 1 



i 2( R 3(fc+l) xR 3(fc+l)) 



Later, the method in Klainerman and Machedon [25] was taken up by Kirkpatrick, 
Schlein, and Staffilani in [23] , where they studied the corresponding problem in 2d, 
and Chen, Pavlovic and Tzirakis [4, 5, 6], in which they considered the Id and 2d 
3-body interaction problem and the general existence theory of hierarchy 1.4. 

We are interested in the case r] l 7^ 0. So we study the Gross-Pitaevskii hierar- 
chy with anisotropic switchable quadratic traps. That is a sequence of functions 
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J <y( fc )(T, yk',yi) f , where ret, Yk,Yk *= ^ nk , which are symmetric, in the sense 
I J fe=i 

that 

7 (fc V,yt;yt)=7 (fe) (T,^;y!) 

and 

7 (fe) (^ I y CT (i) J y<T(2),-,y ( r(fe);y^(i),y^(2),-,ya(fe)) = 7 (fc) ( T ,yi,y2,-,yfc;yi,y 2 ,-,yfc) 

for any permutation a, since we are considering Bosons, and satisfy the anisotropic 
switchable quadratic traps Gross-Pitaevskii infinite hierarchy of equations: 



k 

\h^(t) + \h^(t)) 7 « = &o (V^ 1 )) . (1.6) 

k J 3 = 1 

In the above, Bj^+i = Bjk+i ~ ^jk+i are defined as 

Bl k+l {l {k+1) ){r,yt^ k ) (1.7) 

= J ^(yj -y/c+i)^(yj -yUi)7 (fe+1 H^,y^K;y'fe+i) rf yfe+i d yfc+i 

Bh +1 ( 7 (fe+1) )(r,yt;yt) 
= / /5(y--y fe +i)5(y;-yUi)7 (fe+1) (r,y^l;yC)^yfe+i^ + i- 

These Dirac delta functions in Bj^+i are the reason we consider the collapsing 
estimates like estimate 1.5. 

When the initial data is a BEC / factorized state 

7 «(0,yt;yl) = n^o(y,)^o(y;-), 

3=1 

hierarchy 1.6 admits one solution 

k 

7 (fe) (r,yt;yt) = n^(r,y j )0(^), 

3=1 

which is also a BEC state, provided <f> solves the n — d Gross-Pitaevskii equation 

1 2 
id T (j)~ -H y {r)4> = b \<f>\ 4> 

0(0, y) = o (y). 
Hence we would like to have uniqueness theorems of hierarchy 1.6. 



1.1. Main Auxiliary Theorems. To obtain Theorem 1, we need the auxiliary 
theorems in this subsection which are of independent interest. We show them in 
3d as well. On the one hand, the general idea for the 2d case is derived from the 
higher dimensional case. On the other hand, the 2d and 3d cases are dramatically 
different when they are viewed in the context of Theorem 1. We will explain this 
difference between the 2d and 3d case in Section 7. For the moment, notice that the 
uniqueness theorems in 2d and 3d address two different Gross-Pitaevskii hierarchies 
which stand for the two sides of the lens transform. Also, we currently do not have 
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a 3d version of the 2d convergence / Theorem 1. We state our auxiliary theorems 
regarding different dimensions separately for comparison. 

First, we have the following collapsing estimates which generalizes estimate 1.5. 

Theorem 2. (3*n-d optimal collapsing estimate) Let n — 2 or 3, write 



dxf 



where the L] oc functions a\ satisfy 

ai ^ c > a.e. 
Assume u(t, xi, X2, x 2 ) solves the Schrddinger equation 

iu t + L X1 (t)u + i X2 (t)u ± L x - 2 (t)u = in M 3n+1 
u(0,xi,x 2 ,x' 2 ) = /(xi,x 2 ,x 2 ), 

then 



(1.8) 



|V X | 2 u(t,x,x,x) 



dxdt < C 



|V X1 | 



Theorem 2 is a scale invariant estimate when ai = 1 hence it is optimal. In 
fact, it holds for all n ^ 2. The proof is different for n = 2 and n ^ 3. We name 
the third spatial variables x 2 to match the uniqueness theorems. We point out 
that Kirkpatrick, Schlein and Staffilani proved the almost optimal result for the 2d 
constant coefficient case in [23]. Some other collapsing estimates were attained in 
[7, 19]. 

1.1.1. 2d Auxiliary Theorems. Theorem 2 is the key to show the following unique- 
ness theorem. 

Theorem 3. (Uniqueness of 2d GP with time- dependent coefficients) Let L Xfe be in 

r . — > 00 

Theorem 2 and Bj^+i be defined via formula 1.7. Say < u^ k '{r, y^; y k ) > solves 

I J k=l 

the Gross-Pitaevskii hierarchy with variable coefficients 



idt + L^t) 



V-(t) u 



(fc) 



6o£iW (" (fc+1) ), 



subject to zero initial data and the space-time bound 



-T 



k 

n 

3=1 



v x . 



S j;fe+lU ( fe+1 )(i,-;-) 



dt < C k 



L 2 (R 2k xR 2k ) 



for some C > and all 1 < j < fc. TTierc Vfc, t e [0, T], 



k 

n 

i=i 



0. 



L 2 (R 2fc ; 



In contrast to the standard Elgart-Erdos-Schlein-Yau program, we do not need 
a uniqueness theorem regarding the Gross-Pitaevskii hierarchy with anisotropic 
switchable quadratic traps (hierarchy 1.6) to establish Theorem 1. It is enough to 
have Theorem 3 which has no quadratic potential inside. At a glance, the analysis of 
the above hierarchy based on the Laplacian is unrelated to the hierarchy 1.6 based 
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on a Hermite like operator H y (r). However, Carles' generalized lens transform [3] 
links them together. In fact, the generalized lens transform preserves L 2 critical 
NLS and thus the 2d Gross-Pitaevskii hierarchies. The specific version of the lens 
transform we need is in Section 4. 



1.1.2. 3d Auxiliary Theorems. As mentioned before, the uniqueness theorem here 
addresses a different hierarchy from Theorem 3. Of course we can prove a 3d version 
of Theorem 3. However, the disparity between the 2d and 3d case renders such a 
theorem of little value because the lens transform does not preserve the 3d cubic 
NLS. See Section 7 for detail. 
We consider the norm 



4 fc) 7 (fc) (r,-;-) 



L 2 (R3fcxR 3fc ) 



(1.9) 



in which 



4 fc) = (nJ=i^(r)P y ;.(-r); 



Py(r) 



'd Vl 

a 



^2(T)^+/3 2 (T)y 2 

^3(t)4+/3 3 (t)2/ 3 , 



where /3, solves 



Mr) +Vi(r)/3 l (r) = 0,0,(0) - 1,0,(0) = 0. 



(1.10) 



The operator i/3 ; (r)^- + /3,(r)y; was introduced by Carles in [3]. Lemma 3 and 
relation 5.2 indicate that the norm 1.9 is natural. That is because this operator is 
in fact the evolution of the momentum operator — iV. We will compute it in the 
appendix. 

Through a specific generalized lens transform (Proposition 3) we produce the col- 
lapsing estimate which is the key estimate to our 3d uniqueness theorem regarding 
hierarchy 1.6 when n = 3. 

Theorem 4. Let [s,T] C [0,Xb] and /?, be defined through equation 1.10, assume 
-y( fe + 1 )( T; y fe+i; y^ +i ) satisfies the homogeneous equation 



(1.11) 



Then exists a C > independent of j, k, s, and T s.t. 



Rl k) B 



■j,k+i 



L 2 ([s,T] xR 3fc xR 3fc ) 
-1 

#(fc+l) 7 (fc+l) 



i 2( R 3(fc + l) xR 3(fc+l)) 

where the t on the RHS of the above estimate can be chosen freely in [s,T], 



From Theorem 4, it follows 
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Theorem 5. (Uniqueness of 3d GP with anisotropic switchable quadratic traps) 
Let < ^ k '{T,yk]y k ) \ solve the 3d Gross- Pitaevskii hierarchy with anisotropic 

I J k— 1 

switchable quadratic traps (hierarchy 1.6 when n — 3) subject to zero initial data 
and the space-time bound 

dr€C k 



4 fc) ^, fc+ i7 (fc+1) (r,-;-) 



; xR 3fc ) 



(1.12) 



for some C > and all 1 < j < fc. Then Vfc, r e [0, T ], 



= 0. 



Remark 3. It is currently unknown how to show directly that the limit of 7^ in 
3d satisfies the space-time bound 1.12. 

1.2. Organization of the Paper. We show Theorem 2 for n = 3 first in Section 2. 
Utilizing the same scheme, we prove Theorem 2 for n — 2 in Section 3. Compared to 
[25] which uses the approach in the Klainerman-Machedon null form paper [24] , the 
proofs of Theorem 2 here are closer to Beals and Bezard [2] which is a simplification 
of [24] in the sense that duality takes the place of convolution with surface measures. 

In Section 4, we lay down the tools, a generalized lens transform and its related 
properties, involved in establishing Theorems 4 and 5 whose proofs are in Sections 
5 and 6. Theorem 3 follows from the same procedure. 

In Section 7, we put together the generalized lens transform, Theorem 3, and 
the estimates in Kirkpatrick-Schlcin-Stafhlani [23] to establish Theorem 1. We also 
explain the differences between the 2d and 3d cases there. 

In the appendix, we present an algebraic explanation of the generalized lens 
transform, one of the vital tools in this paper. 

1.3. Acknowledgment. The author's thanks go to Professor Matei Machedon 
and Professor Manoussos G. Grillakis for the discussion related to this work and 
pointing out to him the connection between the generalized lens transform and the 
metaplectic representation, to Professor Remi Carles for sharing the history of the 
lens transform with us, to Mr. Kwan-yuet Ho for telling the author about [9], and 
to Miss Victoria Taroudaki for translating the abstract of the paper into French. 

2. Proof of Theorem 2 when n = 3 / 3*3d Collapsing Estimate 

We will make use of the lemma. 

Lemma 1. [25] Let £ G M 3 and P be a 2d plane or sphere in M 3 with the usual 
induced surface measure dS. 

(1 ) Say < a, b < 2, a + b > 2, then 

dS(rj) . C 



(2) Say e = jr:, then 



I * 1 a 1 \b 

dS(rj) 



a+b-2 ' 





1* i2-e 1 i2-e 


l-v 


\£-v\ \v\ 



c 



Both the constants in the above estimates are independent of P. 
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Proof. See pages 174 - 175 of [25]. 

By duality, to gain Theorem 2 when n = 3, it suffices to prove 



□ 



/ 



Let 



| V x | u(t, x, x, x)h(t, x)dxdt 



ff* ai {s)ds 



|v Xl v X2 v x ,/|| 2 . 




f*a 2 (s)ds 

V Jo a 3( s )dsj 

then it brings the solution of equation 1.8 

Accordingly, the spatial Fourier transform of | V x | u(t, x, x, x) is 
Kil / e < ((«i-€»-€i) T ^(€i-C a -«i)+€i'^C a ±(€i) T ^^)/( €l - & 
which allows us to compute that 

r 

/ |V X | u(t, x, x, x)h(t, x)dxdt 



h(t,t 1 )dtdZ 1 dZ 2 dZ' 2 



(spatial Fourier transform on h) 



J ^ |€ 1 |e < ((^-^- c T ^(€i-«2-€i)+«=^ ± (€0 TA '«i)ft(t,€ 1 )df 



- € 2 - £'2, £2, £'2)d£id£ 2 d£' 2 



^ 1(h) ||V Xl V X2 V x </|| i2 (Cauchy-Schwarz) 
where 

- 12 j- e i((C 1 -« 2 -Ci) T A t (c 1 -« 2 -€i)+€j'A t € 2 ±(€i) T A t «i)^ t)€i)(ft 



/<*)-/- 



|£l ^2 €2 1 I&2I 1^2 1 



"^1^2^2- 



So the target of the rest of this section is to show 

I(h)^C\\hf L2 . 

Noticing that the integral 1(h) is symmetric in |^ 1 — £ 2 — £ 2 | an d 1^2 1 > we deal 
with the region: |£j — £ 2 — £2! > 1^2 1 only. We separate this region into two parts, 
Cases I and II. 

When the "±" in equation 1.8 is " + ", Case I is sufficient. To show the estimate 
for " — " , we need both Cases I and II. 

Away from |£j — £ 2 — £ 2 | > |£ 2 |, there are other restrictions on the integration 
regions in Cases I and II. We state the restrictions in the beginning of both Cases I 
and II. Due to the limited space near "J", we omit the actual region. Please keep 
this in mind during reading. 
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2.1. Case I: 1(h) restricted to the region |£' 2 | < 1^2 1 with integration order 

d£ 2 prior to d£' 2 . Write the phase function of the dt integral inside 1(h) as 

(€i - € 2 - ^ (€1 - €2 - £0 + €^€2 ± (&) T A t & 



+ 2 €2 - 



€2 - 



The change of variable 



leads to 

m . / 



l€il 



(Si-^) T ^t(Si-^) T 



+2Ci'A t « 2 ±(Ci) T >l«Ci)ft( t ^ 1 ) (ft 



^2 2 



l€il 



^2 + ^ 



±U 2 ) T ^2- 

(2.1) 

; 

-^1^2^2 



i(2 



^2 — 



e 1 £1 I2 

£,2 H 2"" 



(Si-i 2 ) T A t/ ( Sl -j 2 ) T 



where 

J(A)(t,€i) = 



l£l 2 e i2&A t S 2e -i2gA t ,t 2 



£ €1 ii 

S2 ~ 



£2 + ^ 



(S 1 -; 2 ) T (A t -A t/ )(S 1 -g 2 )^^ / 



±(&) (^-A t ')&)h(t>,£ 1 )dt l dZ 2 d£ 2 . 



Assume for the moment that 



/ 



J(h)(t,^) dt^C h(;£ 



Li 



2 



with C independent of h or £ 1; then 

Hence we end Case I by this proposition. 
Proposition 1. 

J |J(/)(t,€i)| 2 dt<Cr||/(.,^)||^ 

w/iere C is independent of f or 

Remark 4. To avoid confusing notation in the proof of the proposition, we use 



/(i'>£i) ^° replace 

Proof. Again, by duality, we just need to prove 
J J(f)(t^ l )W)dt 



<C1/(.,£J|| Lf N| L? . 
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For convenience, let 



Then 



J J(f)(t,^)g(t)dt 

|£ |2 e i2^A t € 2e -i2^A t ,€ 2 



<; £1 I2 

S2 — g - 



£2 H — 



— (e- i * (t ' l€l ^ ) /(t',€i)) (e- < *(*.€i.tt) fl (t)) dtdt'dZ 2 dH' 2 



I 



^f e 2i£A t £ 2 (e-W*i,?2)g(tj) dt) (f e" 2 ^'^ (e-W^&fit',^ dt') \^\ 2 d^' 2 



S2 2 



t 1 ii ii 

S2 H 2 



131 





dt 






€2- 


€i-€ 2 
2 


2 


t 1 €1 €2 

S2 + 2 


2 



To deal with the dt and eft' integrals, for every fixed £ 2 > let 



di 2 



then 



du o ai(t)g2.i+fl2(*)g2,2+°3(^2,3 . „ 
* ^ ^/ 



which provides a well-defined inverse t(u). 
Consequently, the integral 



J e 2igA t $ 2 (e-i4>(t4i,i' 2 ) g ^ dt = J e-»l«2l 2 ^ e -**(*(«).€i.Ci) fl (t( u )) 



du 



du, 



is indeed the Fourier transform of 

G(u) = e-WM&MgWu)) 
This is well-defined since 



dt 
du 



f \G{u)\ 2 du = f 

JR JR 



e -i<Kt(u),€i,£2)g(t(u)) 



dt 
du 



du 



= I W)Y 

JR 



dt 
du 
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Hence 



J e 2i&A t £ 2 ^4>(t^ 2 )g(t)^j dt f e - 2i £ A t'^ (e-W^^fi?,^ dt' 



f \Hi\ 2 d& j 

1*2 



lll! 2 <2 



t _ §1 §2 

S2 2 



t Lzga 

S2 2 



2^2 



t i $l~$2 
S2 + 2 



d€ 2 









t Ci-€i 
^2 2 


2 


£ 1 €1 €2 

«2 + 2 



p 2 dpder (spherical coordinate in £ 2 ) 



ll! 2 ^2 
2 — SU P 

P 



l*2| 

(Holder in p) 



V 



p 2 d(T 



< ^ll/(-^l)llL ? Nb<! 

However, 

/ 



t Lz3 

S2 2 — 



l€l 



?2 H 2 — 



/|<V) 









2 
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t €i-€i 

S2 2 




t 1 €1 €2 
S2 + 2 



#2 



p 2 der 



t _ £1 li. 

^2 2 



£2 + ^ 





£2- 


«l-«2 

2 


2 

(icr 


^2 2 






2 |* 2 | 2 



#2 



(Reverse the change of variable in formula 2.1. 

i If!, 



£il / : ^ l2 + 2g su p 



1^2 



V 







Ci-€i 
2 


2 


t €l-«2 
S2 2 




£l *2 *2 


2 - £ |* 2 | 2 - £ 



It' l 2+2e It 


>/ |3-2e 


|S2 SI _ 


S2 



(Second part of Lemma 1) 



^ C. 

In the above calculation, the <x in the first line lives on the unit sphere centered 
at the origin while the cr in the second line is on a unit sphere centered at ^ 1 2 ^ 2 . 
We use the same symbol because Lebesgue measure is translation invariant. 

Thus, 



<C\\f(;^)\\ La Jg\\ Li 



□ 
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Remark 5. Because the integral 1(h) is also symmetric in £ 2 an d €2 when the "± " 
in equation 1.8 is "+", we have acquired the estimate in that case. In Case II, we 
will assume that "± " is "— ". 

2.2. Case II: 1(h) restricted to the region |£ 2 | > |£ 2 | with integration order 

d£' 2 prior to d£ 2 - This time we write the phase function to be 

(€1 - €2 - &f At (€1 - €2 - €2) + €2^2 - {&Y M' 2 
= (€1 - £ 2 f A t - £ 2 ) - 2 («! - £ 2 ) T A t £ 2 + fjA t £ 2 



0(t,€i,€2)-2(€i-€ 2 ) T ^€2- 



and let 

j(£)(t,€l) 



|^ 1 | 2 e -2i(€ 1 -€ 2 ) T A t ^ e 2i(£ 1 -€ 2 ) T A t ,^ 
|£l ~^2 - ^ 2 | 2 l£ 2 | 2 l&l 



e -**(*'.€i.fe) e -i*(*.Ci.€i)ft(f,^ 1 )dt'^d€ 2 



Again, we want to prove 
Proposition 2. 

w/iere C is independent of f or £ 1 . 
Proof. We calculate 



l 2 ? 



2 l> |2 I*/ |2 



|£l ^2 €2 1 l£ 2 l |C: 



/ 
/ 

r \ti\ 2 d4 2 r dj' 2 

J I £9 J £1 — £ o — £0 



(e-*< t '^>/(* , ,*i)) (e-**(*.€i.&) ff (t)) dtdt'd&dt, 



|£l ^2 ^2 1 (£2! 1^2 1 



J e M^ 2 fA t ,^ (e-^'^fit',^)) dt' 



£ 2 r |£i - £ 2 - £ 2 r 1^2, 

Fix £ 1 — £ 2 and £ 2 , write 

e -2 4 (€ 1 -€ 2 ) T A t « 2 ^ e -«0(t,€i,€i) 5 ( f )) dt = J e -m,-^ T A t i' 2 (e-i<Kt,ei,e a ) g (t)) dt 



1 



where uj = (wi , uj 2 , W3) is a unit vector in M 3 . Without loss of generality, we assume 

max{|cji| , \uj 2 \ , |w 3 |} = |wi| 

which implies 

75 * * L 
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Let us further assume that wi > (the proof works exactly the same for the wi < 
case), then we can write 

& = (x,0,0) + (0, yi ,y 2 ) 



u(t) = 2u>\ I a\(s)ds. 
Jo 



Again u is invertible with 



dt V3 



So we have 



/ 



e -2 l (€ 1 -€ 2 ) T A t ^ ^ e -i4>(t,^' 2 )g^ dt 



T A t ,i' 2 (e-iHt^^g^dt 
e -»u(«i|£i-£ 2 k) ^ e -2i|C 1 -« 2 |(O,a;2^3) T A t( „ ) (O^ 1 ^ 2 ) e -i0( t ( u )^ 1 ,^) ff (- i (- M ^ ^_ ^ 



du 



where 



G(u) = e - 2j l«i-€2l(O,W2,W3) T A t( „ ) (O > yi,y2) e -i0(t(«),€i4 2 ) 5 ^( M )) 

which still has the property that 

V3 



dit 



Just as in case 1, this procedure hands us 

J j(m,z 1 Mtjdt 



£l £2 €2! 1^2 



dxdyidy 2 



2 1 w |2 
2 



S$ C 



|€i-€2-€2||€2 

dxdy\dy 2 
111 - £2 - €2 1 1^2! 



sup 



|^ - £ 2 | |^ - £ 2 | z, ^) ) {^d£ 2 

/ l?2l 



2 



/ |^^i) 

J 2 



G(x)F(x, £i) 
dy\dy 2 



d£ 2 



72 "S2 



£l l 2 I2 I2 



G(x) dx I d£ 2 (Holder in x) 



Kil 



I2III21 



/ /" dyidy 2 \ 

N |fc-&-«Sf|tt|'J * 
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The first part of Lemma 1 and the restrictions that — ^ 2 — 1 > II2I an d 

1^2 1 < ll 2 | Sn0W 

^ / 111 I 2 / f dyidy 2 \ , t 



2|€i -€ 2 ii€ 2 r^ J |ii-i 2 -i; 



^ c 



llll'^ 



3-2e i2+2e 



2i€i-€ 2 r" e ie 



which finishes the proposition. □ 

3. Proof of Theorem 2 when n = 2 / 3*2d Collapsing Estimate 

By the proof of the n — 3 case in Section 2, we only need to show these two 
estimates: 

Case I Under the restrictions \£ t - £ 2>old - £' 2 | > || 2if)W | and II2I < | €2 
have 

111 I (f d(T (^2,neJ 



Hi! / d°-(€ 2 ,n e J \At>*n 

/ T77T SU P I / J 77- Nl2 < C 

11 y S»2,netu 2 S»2,new ' 2 / 

where £ 2 neu> and l2.oid arc related by formula 2.1 and we write 

^2,new=P (T with O- e S 1 . 

Case II Under the restrictions — £ 2 — l 2 | > II2I an d |l 2 | > ll 2 |> we have 



H1-I2IH2I v - •/ H1-I2-I2IHL 

where £ 2 = (x,y). 

Lemma 1 plays an important role in giving the corresponding estimates in Section 
2. In the 2d case, the subsequent lemma provides its replacement. 

Lemma 2. Let ^el 2 and L be a Id line or circle in R 2 with the usual induced 
line element dS . 

(1) Say < a, b < 1, a + b > 1, then there exists a C independent of L s.t. 

dS( V ) < C 



L\i-n\ a H b II 



b ^ i^io+6-l " 



(2) Let s = then 



s w U 1 ii-^ii+^-y ^ in 2-26 ' 

Proof. We will show the second part in the end of this section. The first part shares 
the exact same proof with Lemma 2.2 in [25]. □ 
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3.1. Proof of Case I. The change of variable 2.1 turns the restrictions into 

€i-€2 



£2,7 



2 

€i-*' 2 



£l - ^2,oW - ^2 > S^.oid > €i 



2 > 



»2.oZd 



>l£|- 



Notice that ^2,new = P <J -> we m f &c t have 



/It/ 



»2,netu 2 



j 1 €1 £ 2 

*2,nett? 2 



rf<T (^ 2 



; _ §1 ii 

*2,new 2 



1 



1-e 



S2,neui ' 2 



#2 



, . ._ 2 7^2 (Second part of Lemma 2) 

^2 SI — £2 



3.2. Proof of Case II. Recall that £ 2 = (x,y), we estimate 



dy 



ICi I 



l€i-€ 2 ||€: 



il+2e 



sup 



i2-2e i e ii+2e ^2 (First part of Lemma 2) 
J I4i _ ?2I ~ l?2l 



< C. 

3.3. Proof of the Second Part of Lemma 2. Due to 

we can separate the integral as 



sup 



da{r)) 



il-e 



< sup 



upf/ 

77 1 yJS 1 and |£-tj 



+ sup 

M V-'S 1 and |C+r,|^ili 



We will only show 



sup 
M 



d<r(v) 



C 



and |C+r,|^J|i |£ - * ' \$ + ^ |£| 2 & 



since the other part is similar. It is clear that 
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As 



sup 



rfcr(f7) 



and |C+r,|>i|l l^-^ 1 "11 + " , 



c 



7— SUB 

Zt 6 \v\ 



(3.1) 

Rotate S 1 such that | is on the positive x axis, then write 77 = pe %9 for (p cos 9,psin 9) 
and observe: 

• Whence [0,f]U[^,27r], 

\pe ie - (|€|,0)| > |€||sin0| 

because |£||sin0| is the distance between the point (|£|,0) and the line 
{angle — 9). 

• Whence [f,4f], 

1^- (|€|,o)| > |€| 

because pe* 6 — (|£| , 0) is the longest edge in the obtuse triangle which con- 
sists of pe l \ (HI ,0) and pe w - (|£| ,0). 

Insert these two elementary observations into estimate 3.1, we have 



v\ V^ 1 



sup 

I 



and \t+n\> ] -$- ll-^l 1 £ \^ + V\ 1 £ , 



C 

-Tzjsup 



In] KJs 1 IZ-vl 1 ' 



C 



il-e 



C 



sup 

|0 



d0 



|pe*-(|€|,0)| 



l-e 



2 sup 



d0 



p V> |peW-(|€|,0)| 



l-e 



c 



- d6 



|l-e 



d0 



llll sin 1 



1 l-e 



i2-2e ' 



To show the other part, namely 



AS 



sup 

I 



dcr{n) 



and |€-77|>M ll-^l 1 "Il + ^l 1 E ^ 

one just needs to notice 



C 



2-2e ' 



\t + v\ 



(l€l,o)- pc 



:(0+ir) 



then one can proceed as above. Therefore we conclude the proof of the second part 
of Lemma 2. 



4. The Lens Transform / Preparation for Theorem 4 

From now on, we enter the proof of Theorems 4 and 5. We set n = 3 until 
Section 7. In this section, we set up the tools involved in the proof of Theorem 4. 
We build the lens transform we need and state the related properties. For simplicity 
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of notations, we write U^ k+1 \r; s) to be the solution operator of equation 1.11 and 
C/ y (r; s) to be the solution operator of 

1 

2" 







u(s,y) = u s (y). 
i.e. U (k+1 \r; s)jo k+1 ^ solvesequationl.il. By definition, 



UW(T;s) = H(u yj (T; S )U y ,(-T;-s)). 



3 = 1 



To be specific, we need this version of the generalized lens transform: 

Proposition 3. There is an operator T x (t) which satisfies the hypothesis in The- 
orem 2 such that 

C/( fc + 1 )(r;0) 7 [ ) fe+1) 



fc+i 

= n 



n 

u (k+1 \ 



T) (M 2 -k<r)\ 



Mr) 

"lW J/1,1 J/1,2 J/1,3 



/ 



J/fc+1,1 J/fe+1,2 J/fe+1,3 



2/1,1 2/i,2 2/1,3 



Vk+1,1 Vk+1,2 Vk+1,3 



Piir)' /3 2 (tY /! s (t) /WftW/W 

in [—To, To], w/iere a/ and ft are defined as in Claim 1, and u ( fe + 1 ) x^+i ; x' fe+1 ) 
is i/ie solution of 

(id t + L^(t)-L^{t)^u( k+ ^ = OmMW 

^)(0,^;xC) = 7^ +1) . 
The proposition will be a corollary of a sequence of claims. 

Claim 1. Assuming Conditions 1 and 2, for I = 1,2,3, the system 
ai(r)+Vi(T)ai{T) = 0, a,(0) = 0, aj(0) = 1, 

ftW + iiWAM = o,ft(o) = i, ft (o) = o. 

defines an odd ct ; and an even ft G C 2 (R) wii/i //ie following properties 
(1) ft is nonzero in [— T ,T ]; 
(%) T/ie Wronskian of ai and ft is constant 1 i.e. 

d,(r)ft(T)-a / (r)ft(r) = l; 

^ T/ie odd function 



(4.1) 



u;( T ) 



is invertible in [— T ,T ] because 
vi{r) = 



a;(r) 
A(r) 

2 > in [-T , T ]. 
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Proof. We show (1) only since all other statements are fairly trivial. 

Suppose Pi{tq) = for some tq in [—To, To] then /3;(— to) = via f$ x is even. Of 

course r 7^ because /3;(0) = 1. Notice that cos (V-\/sup T |t7 z (t)|^ is a nontrivial 
solution of 

v(t) +sup|?7,(t)|w(t) = 0. 

r 

Since cos ^V-\/sup T |^(t)|^ is not a multiple of cos ^r-\/sup T |?7/(t)|^ must have 

at least one zero in [— t ,to] due to the Sturm-Picone comparison theorem. But 
this creates a contradiction. □ 

Though Claim 1 is elementary, its consequences lying below make our procedure 
well-defined. 

Definition 1. (A reminder of the norm) Let f3 l be defined via equation We 
define 

Py(r)= ^ 2 (r)g|+/3 2 (r)y 2 
and 

R k r = UU P yAr)Py' 3 (-r)- 
Lemma 3. P y (r) commutes with the linear operator 

id T - l -(-A yk+ri {T) |y fe | 2 ). 

Moreover, 

P y {r)U y {r;s)f = U y {r;s)P y {s)f. 

Lemma 4. Say Ki(t,x ,y ) is the Green's function of the Id free Schrddinger 
equation 



2dx 2 



then 



Uy(r;0)uo = I TT " / IT Kl ( oVt' ~aT^' Vol) ) "0(2/01,2/02,2/03)^01^02^03, 



(4.2) 

vaZzd m ifte interval [— T, T] m w/iic/i 77; are Lipschitzian and /?j(t) 7^ 0. 

Proof. Carles computed the isotropic case of formula 4.2 in [3]. We include a proof 
of Lemmas 3 and 4 using the metaplectic representation in the appendix. □ 

We can now prove Proposition 3. On the one hand, via Claim 1, we can invert 
t(r)=v 1 (r) = ^ in [-T ,T ]. 
Therefore, the integral part of formula 4.2 

<H*,x) = J (K 1 (t,x 1 ,y 01 )K 1 (v 2 {v^ 1 (t)),X2,yo2)K 1 (v 3 (v^ 1 (t)),x 3 ,y 03 )) u (y 01 ,y 02 ,y 03 )dy 01 dy 02 dyo3 
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in fact solves 



(id t + Lx(t)) </> 
<X0,x) 



in 



x [-v^(T ),v^(T Q )\ 



where 



L x (t) 



2fe 2 2 / 3 2 ( ur 1 (t))fe 2 2 / 3 2 ( wr 1 (i))feI' 



On the other hand, plugging — t into formula 4.2 yields 

3 

L/ y (-r;0) Uo = (I] 




n Kl (~in~^>~a^I\>yoi) ) u o{yoi,yo2,yo3)dy idy 2dyo3 



i=i 



Whir)- 



because ai and /3 Z are odd while /3 ; are even. 
Whence in [-T ,T 1 



fc+1 



7^ 



fe+i 

n 

i=i 



/ 3 ,.to(M^/)\ 

d g 07m 2 



n 



(fc+l)/ Q; l( T ) 2/1,1 2/1,2 2/1,3 i/fc+1,1 2/fc+l,2 j/fc+1,3 

U { (3i(ry ^(tY f3 2 (ry (3 3 (ry-> ^(r) ' /3 2 (r) ' /3 3 (r) ' 
2/i,2 2/i,3 ffe+1,1 2/fc+i,2 2/fe+i,3N 



if u( fe+1 )(t,x fc+ i;x^ +1 ) solves 



^ t + L^(i)- J L^->(t))^+ 1 ) = OinR^xI-Dr 1 ™,^ 1 ^)] 



(fe+i) 
7o • 



^ +1 )(0,x^t;x' fe+1 ) 
At long last, define 

LJfi), when t G [-ur 1 (7o),ur 1 (7o)] 



ix(t) = { 



Lx^r 1 ^)), when t ^ -L-r 1 ^) or t s$ -^(To) 



then we obtain the desired variant of the generalized lens transform i.e. Proposition 
3. 



5. Proof of Theorem 4 

Without loss of generality, we show Theorem 4 for Bj k+1 in Bj^+i when j is 
taken to be 1. This corresponds to the estimate: 



J* J 

J s JR 



R 3k xR 3fc 



4 fc) 7 (fe+1) (r,yt,yi;y;,yi) dfld y ' k 



(5.1) 



\ re [ ' T ol^2 / ./r3(m-i) x r3(m-i) 
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Vr e [s,T], if 7( fc+1 ) satisfies equation l.ff. 
By Proposition 3, we compute 



4 fc) 7 (fe+1) (T,ylyi;yt,yi) 



(5.2) 



3 iiiiil 



(\^\ 2 MA 2 ) \ 



n 

1=1 



nJ v - Vx ;)J u(t+1)( ?|y' xtxi; ^ Xl 



if we let 



because of the relations 



3,1 



- ^ - ft (r) " 



dy 



2 +ft(r)y^ 2 = 0, 



5 



dyj,i dx jt i ' 



Consequently, 

rT 



[ dr [ i?W 7 (,s+1) (r,yt,yi;yt,yi) 2 dytdyt 

JR 3fc xE 3fc 



fe+1 



dytdyl 



\l=2 

-l 



dr r 

(Pi(r) f k 

inf TT/3f(r)) f°° dt f 



inf TT/3?(r) 



n(Vx J V x ,)j^)(^|,4,x 1 ;2,x 1 

nJv,V x ,)j !1 (^(|||,x1,x 1 ;xtx 1 



dxtdxl 



1 Z=2 



dx£dx fe 



II ( v x, v x;.) I U (fe+1) (t,xt,x i; 2,xi) 



dxtdxl 



where we used the fact that the Wronskian of ai and (i l is constant 1, i.e. 
dt _a 1 {T)P l (T)-a 1 (T)'(3 l ( T ) 1 



dr {Pl{T )f 

as shown in Claim 1. 

A corollary of Theorem 2 tells us that 



Corollary 1. Let L x (t) be the same as in Theorem 2 and vS k+1 ^> verify 
(idt + L^^-L^it)^ =0. 
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Then there is a C > 0, independent of j, k, and vS k+1 ^ s.t. 



< C 



k \ 

II( v ^ v «i)j « (fe+1) (i^,x i; 2,xi; 

fc+i \ 

n( v x, v x;) U (fe+1) (0,^l;x' fe+1 



L 2 (RxR 3fc xR 3fc ) 



i 2( R 3(fc+l) xR 3(fc+l)) 



Whence inequality 5.1 follows. 



6. The Uniqueness of Hierarchy 1.6 

To get Theorem 5, we of course use the Klainerman-Machedon board game 
argument to group the terms. For convenience, we assume b = 1 here. 

Lemma 5. One can express j^(ti, •; •) in the Gross-Pitaevskii hierarchy 1.6 as a 
sum of at most 4™ terms of the form 



L 



or in other words, 



(6.1) 



Herer n+1 = (r 2 ,r 3 , ...,r n+ i), D C [s,Ti] n , fj, m are a set of maps from {2, ...,n+l} 
to {1, satisfying /x m (2) = 1 and f-i m (j) < j for all j, and 

J(z n+1 ,»J = U( 1 \t 1 -t 2 )B 1 ^ 2 \t 2 -tz)B^ ) . 2 ... 

t/ (n) (T n ;T n+1 )B Mm(n+1) , n+1 ( 7 ( n+1 )(r n+1 , •; ■))■ 

Proof. The RHS of formula 6.1 is in fact a Duhamel principle. This lemma follows 
from the proof of Theorem 3.4 in [25] which uses a board game inspired by the 
Feynman graph argument in [14]. One just needs to replace e l (*i-*2)A B by U y (ii; t 2 ), 
ande ! ( il -' 2 > AW by [/« t 2 ). □ 

Let Z? T2 = {(T3, ...,r n +i) | (t2,T3, ...,r n +i) G £>} where D is as in Lemma 5. 
Assuming that we have already verified 



L 2 (R 3 xR 3 ) 



0. 
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applying Lemma 5 to [s,t\] c [0, To], we have 

1 L 2 (R 3 xR 3 ) 



4V / U^(T 1 -T 2 )B 1 . 2 U^(r 2 ;T 3 )B tim{3) , 2 ...dT 2 ...dT n+1 

J D 

£ C/ (1) (t i; t 2 ) ^ i?W J B 1 , 2 C/( 2 )(T 2 ;T 3 )i? Alm( 3), 2 ...dT 3 ...dT„ +1 ^ rfr 2 



L 2 (R 3 xR 3 ) 



(Lemma 3) 



< f 1 / 4 1 2 )S 1 , 2 C/( 2 )(T 2 ;r 3 )i? AJm(3) . 2 ... ( iT 3 ...dT„ +1 
^ / R^B^ 2 U {2 \T 2 -T Z )B^ (m ... 

J [s,Tl]" 

< (ri - s) 



dr 2 



L 2 (R 3 xR 3 ) 

4 1 2 )B i,2^ (2) (T 2 ;r 3 )i? Alm(3)i2 .. 



L 2 (R 3 xR 3 ) 

dT 2 dT 3 ...dT n+1 



dr 3 ...dT n+1 



M m (3),2- 



L 2 (R 6 x 



^ C(n-*)' / 4 2 2 )C/( 2 )(T 2 ;r 3 )i?, 

•'[s.ti]'- 1 

(Same procedure n — 2 times) 

J s 

^ C (C (n - s))^ . 

Let (ri — s) be sufficiently small, and n — > oo, we infer that 



I 2 (T2e[s,Tl]xl 3 Xl 3 ) 

dr 3 ...c?T„ + i (Theorem 4) 



L 2 (R 3 "xR 3 ™) 



ra+l 



4V7 (1) (n,-) 



Similar arguments show that 
we have attained Theorem 5. 



L 2 (R 3 xR 3 ) 

i#°7 (fc) 0v) 



in [s,n]. 

= 0, Vfc,r e [0,T ]. Hence 



L 2 (R 3 xR 3 ) 



7. Derivation of the 2d Cubic NLS with Anisotropic Switchable 
Quadratic Traps / Proof of Theorem 1 

For a more comprehensible presentation, let us suppose 



H y (r) 



d 2 



dy 



3,1 



+ Vi(r)yli 



is the ordinary Hermite operator 



H y = -A y + |y| 2 

in this section to make formulas shorter and more explicit. We will add two remarks 
in the proof to address the small modifications needed for the general case. 

We start by reviewing the standard Elgart-Erdos-Schlein-Yau program in this 
setting. 
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Step A. Observe that, by definition, {7^'} s °1 vgs the quadratic trap Bogoliubov- 
Born-Green-Kirkwood-Yvon (BBGKY) hierarchy 

(*-i(-^ + W*) + 5(-^ + ^0)^ <"> 



N — k , /■ 

+ N J d Vk+i[(VN(yi-yk+i)-V N (y' i -yk+i)) 



3=1 ' 



(fe+l)/ — > / M 

where Vjv(x) = N n ^V (iV^x). It converges (at least formally) to the qua- 
dratic trap Gross-Pitaevskii infinite hierarchy 

(*-5(-^ + W") + K-^ + l^0) 7W (7 ' 2) 



Prove rigorously that the sequence {7^} ^ s com P ac t with respect to the 

weak* topology on the trace class operators and every limit point {7^} 
satisfies hierarchy 7.2. 
Step B. Utilize a suitable uniqueness theorem of hierarchy 7.2 to conclude that 



-> k 



l (k) (r,yt;y[) = l[ ( t ) (T,y j ) < f>(T,y' j ), 
3=1 

where <j) solves the 2d quadratic trap cubic NLS 

^ T 0=i(-A + |y| 2 )0 + 6o0H 2 . 

So the compact sequence {7^} nas on ^y one limit point, i.e. 

k 

7^ fe) ^n^ r '^)^yp 

3=1 

in the weak* topology. Since 7^ is an orthogonal projection, the conver- 
gence in the weak* topology is equivalent to the convergence in the trace 
norm topology. 

We modify this procedure to show Theorem 1. We remark that the main addi- 
tional tool is the lens transform. When H y (r) is the Hermite operator, ai = sinr, 
01 = cost and T < | i.e. the lens transform and its inverse reads as follow. 
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Definition 2. We define the lens transform Ti : L 2 {d^dx k ) — > L 2 (dytdy\.) and 
its inverse by 

f^Wy!;^) = r J- r ^)(tanr,^;i)e-^ ( N 2 -|^r ) 
V / (cost)"* cost cost 

(t-^A it xl-xfl - 1 -yWfarctan/ — S ■ ** = y^ig?) (N 2 -|**| 

l< 7 J (< ' Xfc ' Xfe) - (i + i2) ^ 7 (arctan *' vr+F' 

7] is unitary by definition and the variables are related by 
t = arctani, y fe = Xfc and yj, = 



Remark 6. For i/ie general anisotropic case, we still need the 2d version of Propo- 
sition 3. 

Let us write 

(ij-y fc) ) (t, xfc 2) = 7 (fc) (r, yt; yt) eiIT7 ^ ( ' Xfc ' J^ 1 - := 7 (fc) (r, yt ; yt)^ fc) (*, 4; 2), 
v 7 (l + t 2 )~ 

then we have a more explicit version of Proposition 3. 
Proposition 4. 

(idt + 1 -A^ k - 1 -A^j (V 7 ( fe )) (t, x|; xt) 

- \ (-Ayj + |yt| 2 ) + ^ (-A^ + |yt| 2 )) 7 (fe) (r,yt;yt) 



/i (fc) 



l + t 2 

Proof. This is a direct computation. □ 



Via this proposition, we understand how the lens transform acts on hierarchies 
7.1 and 7.2. 

Lemma 6. (Gross-Pitaevskii hierarchy under the lens transform) {7^} solves the 
quadratic trap Gross-Pitaevskii hierarchy 7.2 if and only if {u^ — Tf 1 ^^ solves 
the infinite hierarchy 

U + - \ A x -,) „« = ^£b £ B, k+l . (7.3) 

.7=1 



In particular, when n = 2, the lens transform preserves the Gross-Pitaevskii hier- 
archy. 
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Lemma 7. (BBGKY hierarchy under the lens transform) {7^} s °l ves the qua- 

dratic trap BBGKY hierarchy 7.1 if and only if < u N = T l ^y N > solves the hier- 
archy 



1 1 

Nl + t 2 



k 



N 



(7.4) 



+ V 1 + t2 E/^ +1 [(M 



) « 



(fc) 

N 



Vn{ ^ttW } 



u { n +1) (*, x fe+ i; x£, x fe+ i)], 



We can now prove Theorem 1. 
7.1. Proof of Theorem 1. 

Step 1. Let n = 2, consider = ^ _1 7^ j which solves hierarchy 7.4. 

Step 2. Write 

V(x) = -J_y(_^=), 



then 



(1 + T 2 ) 1 ^ 



when T < 00 and p > 1. 



Therefore we can employ the proof in Kirkpatrick-Schlein-Staffilani [23] 

to show that the sequence jw^j is compact with respect to the weak* 

topology on the trace class operators and every limit point {it^} satisfies 
the Gross-Pitaevskii hierarchy 7.3. Moreover, based on a fixed time trace 
theorem argument as in [23], for a < 1, we have 



dt 



n(<vx,r<v x ,niw 



€ c k . 



L 2 (R 2k xR 2k ) 



for every limit point {u^ fe H. To be more precise, the proof in [23] involves 
a smooth approximation. We omit this detail here. 

Remark 7. The auxiliary Hamiltonian 

1 N 1 

H N {t) = - £ L Xj (t) + -J2 N 2 "V(N^ (x, - x,)). 



i<3 



which corresponds to the anisotropic quadratic potential case does not lead to the 
conservation of the quantity 
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On the other hand, the following estimate controls the energy. 



N 



+ 



+ ... 



iP N , (^iv(i)) 
CkU N ,{H N {t)) k i, 



Jt ,[H N (t) 



N 



since a\ and a-i, the coefficients of Lx, are C 1 in the context of Theorem 1. Thus 
Gronwall 's inequality takes care of the problem for us as long as we are considering 
finite time. 

Step 3. By Theorem 3 (2d uniqueness) or Theorem 7.1 in [23], we deduce that 

k 

where </> solves the 2d cubic NLS 

1 



id t (f> = --A(f> + b <i) 



Hence the compact sequence j has only one limit point, so 



in the weak* topology. Since is an orthogonal projection, the conver- 
gence in the weak* topology is equivalent to the convergence in the trace 
norm topology. 

Remark 8. It is necessary to use Theorem 3 in this paper for the general anisotropic 
quadratic traps case. 

Step 4. Let (j> solve the 2d quadratic trap cubic NLS 

z9 T 0=i(-A + |y| 2 )0 + 6 o 0|^| 2 , 
then the lens transform of vP^ is 

k 

7 (fe) (r,yt;yt) = n^^)^I)' 

due to the fact that the lens transform preserves mass critical NLS, which 
is the cubic NLS in 2d. 
Step 5. The convergence 

in the trace norm indicates the convergence in the Hilbcrt-Schmidt norm. 
But the lens transform 



Ti : L 2 (d^dx) -> L 2 {d~fdy) 
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is unitary (so preserves the norm) and thus 
Thence we conclude that 7^ converges to 

k 

in the Hilbert-Schmidt norm, which is Theorem 1. 

7.2. Comments about the 3d case. It is natural to wonder what we can say 
about the 3d case using the above method. Visiting Lemma 6 again yields the 
hierarchy 

U + iA^ - Ia^) «« = (1 +i 2 )H £iW (u^) . (7.5) 

Due to the factor (l + t 2 ) 2 , it is difficult to see of what use a 3d version of Theorem 
3 might be. We can certainly give a uniqueness theorem regarding hierarchy 7.5 
with the techniques in this paper. But it is unknown how to verify the space-time 
bound when n = 3 as stated earlier, 

Another possibility to attack the 3d case is the standard Elgart-Erdos-Schlcin- 
Yau procedure, but we presently know very little about the analysis of the Hermite 
like operator H y {r). 

Finally, we remark that it is not clear whether the Feynman diagrams argument, 
the key to the uniqueness theorem in [14] on which [13, 14, 15, 16, 17] are based, 
leads to a 3d uniqueness theorem of hierarchy 1.6 or 7.5, which represent the two 
sides of the lens transform. 

8. Conclusion 

In this paper, we have derived rigorously the 2d cubic NLS with anisotropic 
switchable quadratic traps through a modified Elgart-Erdos-Schlein-Yau procedure. 
We have attained partial results in 3d as well. Unfortunately, when n — 3, we still 
have unsolved problems as stated in Section 7.2. 

9. Appendix: the Generalized Lens Transform and the Metaplectic 

Representation 

In this appendix, we prove Lemmas 3 and 4 via the metaplectic representation. 
The 3d anisotropic case drops out once we show the Id case. Before we delve into 
the proof, we remark that we currently do not have an explanation away from direct 
computations for Proposition 4 or for the fact that the generalized lens transform 
preserves L 2 critical NLS. The group theory proof presented in this appendix only 
shows the linear case: Lemmas 3 and 4. 

Through out this appendix, we consider the metaplectic representation 

pi : Sp (2,R) -> Unitary Operators on L 2 (K). 
which has the property: 

*(U, ;))-'(-^>£)- 
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For more information regarding /i and dfx, we refer the readers to Folland's mono- 
graph [18]. We comment that \i is not a well-defined group homomorphism on all 
of 5p(2,R) , but the fact that it is well-defined in a neighborhood of the identity 
of Sp (2, K) is good enough for our purpose here. 

9.1. Proof of Lemma 4 / the Generalized Lens Transform. 

Proposition 5. Define a and ft through the system 

a{r) + T]{r)a(T) = 0, a(0) = 0, d(0) = 1, 
Ht) + v{t)0{t) = 0,0(0) = 1,0(0) =0, 

and let 

0(t) -a( T y 



B{T) " \-P{t) d(r) 
Assume (3 is nonzero in some time interval [0, T], then \x (B (r)) / solves the Schrddinger 
equation with switchable quadratic trap: 

id T u = ^-idl + r](T)^pjuinRx[0,T] (9.1) 

u(0,y) = f(y)eL 2 (R). 

Proof. We calculate 

d T \ T ^(B(T Q +r))f = (d T \r=0»(B(T +T)))f 

= (d T \ T= oV (B (r + t) B- 1 (t ) B (t ))) / 
= (d T | T=0 M (B (to + t) B- 1 (t ))) /i (.B (t )) / 
= d M (B'(ro)£-Vo))M0B(To))/. 



where 

B'(to)B _1 (t ) 



0(t o ) -d(r )\ /d(T ) a(r ) 

-0(t o ) a(T ) y ^(to) 0(t„) 

f 0(r o ) -a(ro) \ A*^o) a(r )\ 

W^)0(r o ) -r?(T )a(T )y ^0(t o ) 0(to)J 



/ 0(r o )a(ro) - d(T O )0(T O ) N 

^(t ) (d(To)^(To) - ^(To)a(To)) y 

Notice that the Wronskian of a and is constant 1 i.e. 

c*(t)0(t) - <*(t)0(t) = 1. 

So 

d^B'i^B- 1 (t )) = d/x ° 



In other words, 



V(r ) 
= -^(-^ 2 + ^(ro) y 2 ). 



9 r (/x (B (t)) /) = - I (-d 2 + r7(T) y 2 ) („ (B (t)) /) . 

Before we end the proof, we remark that ^ is required for the metaplectic 
representation to be well-defined. □ 
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Through the LDU decomposition of the matrix B, we derive the generalized lens 
transform. The LDU decomposition of the matrix B is 

B(t) = 






Hence we have 

n{B{r))f = n\\ kT) ;||M(r V n" _L) ) H ( C W> I I /> ( 9 - 2 ) 
where 

e iP M^f(y) by (4.25) in [18] 

r/(flfr) by (4.24) in [18] 




•o(t) 9£ 

= e 4 — — / by (4.54) in [18]. 

Due to the definition of fj,, equality 9.2 in fact holds up to a " ± " sign which 
depends on the time interval. However, the LHS and the RHS of equality 9.2 agree 
for sufficiently small r. By continuity, they must agree on the time interval [0,T] 
where /3 0. So we conclude the following lemma concerning the generalized lens 
transform. 

Lemma 8. [3] Assume (3 is nonzero in the time interval [0,T], then the solution 
of the Schrddinger equation with switchable quadratic trap (equation 9.1) in [0, T] 
is given by 

e'^> 2 a r y 

(/3(t)) 5 ^( t ) ^( t ) 
if v(t,x) solves the free Schrdrdinger equation 

id t v = --AiBR 1+1 
2 

w(0,x) = f(x)eL 2 (R). 

The anisotropic case, Lemma 4, follows from the above lemma. 

9.2. Proof of Lemma 3 / Evolution of Momentum. Using the metaplectic 
representation, we can also compute the evolution of momentum and position. 

Lemma 9. The evolution of momentum and position is given by 

P(t) = A1 ( J B(r))oHa,)o( A1 (i?(r)))- 1 =^/3(r)9 y -^(r) 2/ 
Y(r) = ^(BiT^oyo^iBir)))- 1 = ia(r)d y + a(r)y . 
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Proof. Let us only compute the momentum, position can be obtained similarly. 
MB(r))(-^)0*(B(r)))- 1 = H (B(r)) (l 0) ( M (B(r)))- 1 

= (1 0) (B(r)) T (Theorem 2.15 in [18]) 

= (i o) ( -ft T A (~ id A 
[ j \-<*{t) d(r) A y ; 

= -0{r)d v - $(r)y 

□ 

Remark 9. VFe select — «9 y to &e f/ie momentum to match the canonical commu- 
tation relations in Folland [18] which is 

[-id y ,y] = -il. 

The above lemma reproduces the following result in Carles [3] . 
Lemma 10. [3] The operators P(t) and Y(t) commute with the linear operator 

Moreover, 

P(t)U(t;s) = U(t;s)P(s) 
Y(t)U(t;s) = U(t;s)Y(s) 
if we let U v (t; s) be the solution operator of 

id T u = ^d 2 y +f ] (r)^ju l nR 1+1 

u(s,y) = u s (y) g L 2 (R), 

or in other words 

U y (r;s)=n(B(r))n(B(s))- 1 . 
Thence we have shown Lemma 3. 
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